Confinement-induced collapse of a dipolar Bose-Einstein condensate 
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We report on the observation of the confinement-induced collapse dynamics of a dipolar Bose- 
Einstein condensate (dBEC) in a one-dimensional optical lattice. We show that for a fixed interaction 
strength the collapse can be initiated in-trap by lowering the lattice depth below a critical value. 
Moreover, a stable dBEC in the lattice may become unstable during the time-of-flight dynamics 
upon release, due to the combined effect of the anisotropy of the dipolar interactions and inter-site 
coherence in the lattice. 

PACS numbers: 03.75.-b, 67.85.-d 



The stability of many-body systems, either classical 
or quantum, is typically determined by inter-particle in- 
teractions. Crucially, a change in the interaction energy 
balance may drive the system into instability, followed by 
a collapse, as spectacularly exemplified by stellar super- 
nova explosions. Such interaction-induced instabilities 
have been observed as well in degenerate quantum gases, 
where a change in inter-atomic interactions by means of 
Feshbach resonances has been shown to induce the so- 
called Bose-Nova explosion [TJ [2] . 

Recent experiments are exploring the physics of sys- 
tems with significant dipole-dipole interactions (DDI), 
including highly magnetic atoms as chromium [3j |4], 
dysprosium [5^ and erbium [6 , and ultra-cold polar 
molecules [?]• Even for the more complex interaction 
landscape provided by the DDI, the instability and col- 
lapse of a dBEC may be driven by a change in the in- 
teratomic interactions, the remarkable difference being 
that the anisotropy of the DDI leads to different collapse 
dynamics depending on the initial trap shape [HI |9] . 

In this work we show that the collapse of a dBEC can 
be induced as well by a change in the external trap- 
ping potential, while keeping the interaction strength 
constant. Such confinement-induced collapse, contrary 
to the interaction-induced case, relies on the geometry- 
dependent stability of a dBEC [TOl |TT]. Thus, it is a 
general feature of dipolar systems, which cannot be ob- 
served in non-dipolar BECs, nor to our knowledge, in any 
other many-body system. 

We investigate the collapse dynamics of a dBEC with 
a fixed short-range interaction strength, trapped in a ID 
optical lattice. Starting with an initially stable dBEC, we 
drive the system into instability by reducing the depth 
of the lattice potential below the stability threshold that 
we have recently mapped [TP . We show that this change 
in the external confinement induces an in-trap collapse, 
revealed by strong in-trap atom losses. Moreover, we 
show that a stable dBEC in the lattice may become un- 
stable and collapse in time-of-flight (TOP) upon release. 
This TOF-induced collapse is a general feature of dipolar 
BECs in optical lattices, which, as we show, results from 



both the anisotropic DDI and the coherence of the BEC 
in the lattice. As TOP imaging is a major measurement 
technique, especially in lattice gases, the TOP-induced 
collapse demonstrated here is a key issue to be consid- 
ered in experiments with polar lattice gases. 

Our experimental procedure is as follows. We first pre- 
pare a quantum gas of bosonic ^^Cr atoms in a stable con- 
figuration, following the sequence presented in Ref. [TT]. 
The BEC, containing typically 15,000 atoms, is confined 
in the combined potential produced by a crossed optical 
dipole trap (ODT), characterized by harmonic frequen- 
cies I'x.y.z = (540, 270, 470) Hz, and a ID optical lattice, 
with inter-site spacing diat = 534 nm, oriented along z. 
The atomic cloud is polarized by a strong magnetic field, 
also oriented along in the vicinity of a Peshbach reso- 
nance. We make use of this resonance to tune the 5-wave 
scattering length, characterizing the short-range contact 
interaction, down to a = (2 ± 2) ao, with ao the Bohr ra- 
dius. The lattice depth, initially equal to U\nit = 12.6 
(where = h^n'^ / {2md^^^) is the recoil energy, with m 
the atomic mass) is then ramped down to its final value 
U in 100 /is, while keeping the scattering length constant. 
U can be chosen arbitrarily above or under the stability 
threshold, the specific values that we are using in this 
paper being shown in Pig. [l] We hold the system at this 
final configuration for an adjustable time thoid and fi- 
nally switch off all optical trapping potentials to perform 
an 8 ms TOP before taking an absorption image [12 . 

The time-evolution of the system is shown in Pig. |2] for 
increasing holding time for different final lattice depths 
U < t/init- Each snapshot results from the average of 
five absorption images taken after TOP under the same 
experimental conditions. Before averaging, the broad 
isotropic thermal cloud present on each single image was 
fitted by a Gaussian and subtracted from the image. Fig- 
ure |2] therefore shows the atomic patterns resulting from 
the interference of the remaining coherent atoms. 

The system is either stable or unstable in-trap (i.e. be- 
fore release) depending whether U is chosen above or un- 
der the stability threshold, located at around 7 Er for our 
parameters (see Pig. [T]). To distinguish between the dif- 
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Figure 1: Stability diagram of the ^^Cr BEC trapped in the 
ODT and the ID lattice. The red solid line corresponds to 
the stability threshold calculated for the experimental param- 
eters nX. The system is initially prepared in a stable con- 
figuration ([/init = 12.6 ^R, a = 2 zb 2ao), denoted by the 
red square. The blue/dark gray (red/light gray) crosses cor- 
respond to the different values of the final lattice depth U 
chosen in the stable (unstable) region. The inset depicts the 
time sequence of the experiment. 



ferent in-trap dynamics expected for these two regimes, 
we determine as a function of thoid the evolution of the 
remnant fraction^ defined as the number of remaining co- 
herent atoms (integrated over the images in Fig. [2]) nor- 
malized to the total atom number before subtracting the 
thermal cloud (see inset of Fig. [3|. From an exponential 
fit of the time evolution of the remnant fraction, we ex- 
tract the atom loss rate, which serves as an observable 
for the in-trap dynamics. Loss rates are shown in Fig. |3] 
as a function of the final lattice depth U before the TOF. 

We observe a clearly different dynamics depending on 
the value of U compared to the stability threshold. For 
U < 7^R, the atomic cloud experiences an in-trap dy- 
namics characterized by strong atom losses, showing that 
indeed the dBEC becomes unstable in-trap, before re- 
lease. This collapse relies on the anisotropy of the DDI 
and the trap geometry after ramping down the lattice 
depth. 

On the contrary, for U > 7 Eji the system presents 
almost no evolution with the in-trap holding time thoid, 
as it can be seen from the very low loss rates. In addi- 
tion, the sudden release from the lattice after the hold- 
ing time results for all thoid in the usual interference 
pattern formed by a central peak, corresponding to the 
zero- momentum component, and two side peaks, asso- 
ciated with the lattice recoil momentum 2hk\g,t (with 



^lat = 7r/(iiat) [13^. However, in contrast with typical 
interference patterns obtained from non-dipolar BECs, 
the central peak exhibits in our case a clear d-wave sym- 
metry, similar to the one observed in interaction-induced 
collapse experiments [8\ From the absence of evolution 
with thoid and the observation of the d-wave shape of the 
central peak, we deduce that the collapse of the dBEC 
happens during the TOF: the system being stable before 
release, the collapse is therefore induced by the TOF it- 
self [14 . As discussed below, this TOF- induced collapse 
is specific to dipolar gases in optical lattices and relies on 
the non-trivial interplay between the anisotropy of the 
dipolar interaction and the coherence of the system in 
the lattice. 

To examine more closely the condensate dynamics, 
we perform numerical simulations based on the nonlocal 
nonlinear Schrodinger equation: 



+N J dr'l^„t(r-r')|*(r',t)| = 



^{r,t) (1) 



where ^(r, t) is the condensate wavefunction, m the 
atomic mass, and N the initial number of atoms. 
The external potential 14xt(r) = U sin^(7rz/(iiat) + 
m ^.^^ ^ ^(27ri/i)^rf /2 results from the combination of 
the ID optical lattice and the 3D harmonic confinement 
given by the ODT. The interaction energy is given by 
both contact and DDI potentials Fint(r) = ^^^S{r) + 
^1 — 3^^, where the dipoles are polarized along z. 
The non-unitary term proportional to the loss rate L3 
models non-resonant three-body losses, being essential 
for a realistic simulation of the collapse dynamics [2^ [8] . 
L3 is taken equal to 2 x 10~^^ m^/s [8]. The confinement- 
induced instability is then simulated by means of real- 
time evolution of Eq. ([T]), taking into account the whole 
experimental procedure described above, including the 
TOF expansion. 

Results of our numerical simulations are presented in 
Figure [4j where we show snapshots of the collapse dynam- 
ics of the system and the time evolution of the number 
of atoms, both for different final lattice depth U. Here 
we set the origin of the time axis at the end of the lat- 
tice ramp and we let the atomic cloud evolve for 0.6 ms 
in-trap before releasing it. 

We first focus on the extreme case U = OEr below 
the stability threshold. In this case, we observe that the 
atomic cloud shrinks and undergoes strong atom losses 
while it is still trapped. Once released, however, the 
atomic cloud does not suffer any atom loss anymore. This 
shows indeed that the instability of the atomic cloud is 
initiated in-trap and that the subsequent in-trap collapse 
dynamics is associated with strong atom losses. After col- 
lapsing in-trap, the atomic cloud acquires a shape similar 
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Figure 2: Collapse dynamics: series of images (taken after an 8 ms TOF) of the collapsing system for different in-trap holding 
times thoid, at different final lattice depth U above and under the stability threshold (located around 7^r). Each image is 
obtained by averaging 5 absorption pictures, with the thermal cloud removed. 
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Figure 3: Evolution of the loss rate of the remaining coher- 
ent atoms for different final lattice depths U. The loss rate 
strongly increases when crossing the stability threshold. The 
inset shows the time-evolution of the remnant fraction (see 
text for definition) for different final lattice depth [/, from 
which the loss rates are extracted. 



to the one observed experimentally and clearly visible af- 
ter release. The dynamics of the system is different when 
U gets closer to the stability threshold. This is particu- 
larly clear for U = 6.3 £^r: in this case, our simulations 
show that the atom losses start while the system is still 
trapped but only stops well after the beginning of the 
TOF. Therefore, even though the system starts to evolve 
in-trap, it only collapses with a d-wave symmetry during 
the TOF due to the short in-trap holding time. 



In contrast, in the case U > 7£^r, we observe that the 
system is basically not evolving in-trap and undergoes 
only very low in-trap atom losses. This behavior reflects 
the stability of the atomic cloud before the TOF. In ad- 
dition, our simulations show that the evolution of the 
system after release from the trap is a two-step process: 
first the high-momentum components of the wave func- 
tion separate from its zero- momentum component. Then 
the latter shrinks transversally and collapses with a d- 
wave symmetry. This collapse of the central cloud can 
also be observed from the atom number, which exhibits 
a sudden and localized decrease during the TOF (occur- 
ring at around t = 1.45ms for U — 12.6 ^r). Hence, our 
numerical simulations confirm the existence of the two 
types of collapse: an in-trap collapse and a TOF-induced 
collapse. 

We now consider the particular case of the TOF- 
induced collapse. To get more insight into this peculiar 
phenomenon, let us consider a deep enough lattice, such 
that the single-band approximation holds. We consider 
also for simplicity a decoupling of the radial (x, y) and 
axial {z) coordinates. Assuming inter-site coherence, the 
in-trap momentum distribution along the lattice direc- 
tion is given by a series of narrow peaks at ki = 2k\g,tl 
(with / G Z), whose form is given by F (kz), the Fourier 
transform of the spatial envelope F{z) associated to the 
overall harmonic confinement along z. Since the recoil 
energy is much larger than any other energy scale in the 
problem, at the early stages of the TOF (ttof ~ ^/^r) 
particles belonging to different momentum peaks quickly 
move away from each other with a large relative velocity 
On the contrary the form of the wave packet 
associated to each individual peak evolves much more 
slowly, being given during these early TOF stages by the 
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Figure 4: Real-time simulations of the experimental collapse sequence shown Figure [l] The first (second) row shows snapshots 
of the time-evolution of the system with a final lattice depth oi U — 12.6 {U — OEr). Note that the snapshots in the two 
rows are taken at different times, t = marks the end of the lattice ramp and at t=0.6ms, the trapping potential is suddenly 
switched off. The lower panel wall shows the time-evolution of the atom number for different final lattice depths U, above as 
well as under the stability threshold. 



original envelope F{z). The populations of the different 
wave packets are given by A^|VK^p, where Wi = W{ki) is 
the Fourier transform of the on-site Wannier wave func- 
tion. 

Each wave packet evolves then independently. Note 
that, crucially, interactions remain relevant for ttof > 
h/E^i due to the slow expansion dynamics of each in- 
dividual wave packet. As a result, the wave packets may 
either expand indefinitely or collapse, depending on their 
population and geometry (due to the anisotropy of the 
DDI). For our typical parameters, a variational analysis 
following a similar Gaussian ansatz as that introduced 
in Ref. [15 shows, in agreement with our experimental 
results, that the central peak collapses whereas the side 
peaks do not. This results from the different relative 
populations of the momentum peaks, which can be easily 
calculated by approximating the Wannier function by a 
Gaussian. For U = 12£^r, thej"elative population of the 
zero-momentum component \Wo\'^ = 0.61 is much higher 
than th^ relative population of each 2 fi/ciat- momentum 
peak = 0.19, leading to the collapse of the central 

peak only. 

Note that the role of inter-site coherence is crucial in 
the above discussion. Indeed, in absence of coherence, 
the in-trap momentum distribution presents no individ- 
ual peaks, but rather a broad Gaussian-like distribution 
W{k). As a result, after a time scale ~ h/Eji^ the in- 
coherent sum of the expanded wave functions of each 



lattice site results in a single broad wave packet with a 
rapidly growing width. Thus, within this time scale, the 
atomic density drops dramatically, interactions become 
irrelevant and the atomic cloud expands freely. In addi- 
tion, the collapse discussed above crucially depends on 
the anisotropy of the DDI, and the change of the cloud 
geometry in TOF. In particular, non-dipolar BECs, even 
with a < 0, do not collapse in TOF if they were stable 
in-trap. 

Finally, we note that confinement-induced collapse, ei- 
ther in-trap or TOF-induced, can be generalized to other 
dipolar systems as long as the scattering length of the 
system is chosen below add, where add = m iiq ii^ /127i fi? 
is the length scale associated to the DDI, with ji the mag- 
netic dipole moment. In the case of chromium, /j, = 6 fiB 
leading to add — 15 ao and as expected, we do not observe 
any collapse arising when the same experiment is per- 
formed on a dBEC with a scattering length above 15 ao- 

In conclusion, we have shown that a dBEC may col- 
lapse under a change in its trapping potential. Such 
confinement-induced collapse, performed at constant in- 
teraction strength, is in strong contrast with previously 
studied interaction-induced collapses, being a character- 
istic feature of dBECs. We have furthermore shown that 
a dBEC initially stabilized in a ID optical lattice can 
exhibit two different types of confinement-induced col- 
lapse, depending on the confinement configuration before 
release. On the one side, a dBEC may be destabilized 
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while still trapped (in-trap collapse). On the other side, 
a stable dBEC in the lattice may become unstable and 
collapse in TOF (TOF-induced cohapse). 

TOF-induced collapse is a characteristic feature of 
dBECs in lattices of moderated depths resulting from 
the anisotropy of the DDI and inter-site coherence. It 
shows that, contrary to the typical assumption that in- 
teractions do not play any role in TOF or only introduce 
a distortion in the momentum distribution, the TOF dy- 
namics of a dBEC is more complex than expected. Since 
TOF imaging is a basic tool in ultra-cold gases, especially 
in lattice experiments, we stress that the TOF-induced 
collapse demonstrated here will have important conse- 
quences for future experiments on polar lattice gases, 
constituted either by atoms [5l |6] or molecules [7] . 
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